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ON  JETS  PRODUCED  BY  CONICAL  NOZZLES 
by  Ignace  I.  Kolodner 

When  a  steady  stream  of  fluid  is  ejected  from  a  conical 
nozzle,  a  jet  formshaving  the  shape  of  a  hollow  tulip  at- 
tached to  the  orifice  of  the  nozzle.  It  is  also  observed 
that  the  jet  will  oscillate  about  its  equilibrium  config- 
uration, unless  it  is  vented  by,  say,  piercing  it  at  some 
point. 

These  observations  can  be  explained  qualitatively  as 
follows.   The  tulip  shape  of  the  jet  results  from  the  action 
of  surface  tension  which  tends  to  contract  the  jet  toi^rard 
the  axis  of  the  nozzle.   Thus,  for  a  given  set  of  data  there 
forms  an  equilibrium  jet  envelope.   If  for  some  reason  the 
jet  envelope  is  made  to  deviate  initially  from  its  equili- 
brium shape  and  if  it  is  closed  (I.e.,  unvented),  oscillations 
will  be  induced.   Indeed,  a  contraction  of  the  jet  com- 
presses the  air  enclosed  by  it  and  thus  produces  an  out- 
ward force  which  would  reexpand  the  jet  and  cause  it  to 
overshoot  its  equilibrium  position,  after  which  it  would 
be  stopped  and  pushed  in  again.   If  the  jet  is  vented  no 
oscillations  will  occur,  since  the  inside  air,  being  kept 
at  fixed  pressure,  would  not  supply  an  outxijard  force.   The 
main  point  to  be  stressed  here  is  that  surface  tension  alone 
could  not  account  for  the  appearance  of  oscillations.   Of 
course,  a  vented  contracted  jet  vjill  resume  its  equilibrium 
shape  because  the  new  supply  of  liquid  from  the  nozzle  will 
essentially  follow  the  equilibrium  path,  while  the  fluid  in 
the  contracted  portion  will  flow  out. 

The  object  of  this  paper  is  the  quantitative  study  of 
these  phenomena  carried  out  under  the  assumption  of  the 
thin  jet  theory.   In  section  P  the  mathematical  problem  is 
formulated  and  equations  of  motion  of  the  fluid  ejected  from 
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the  nozzle  are  derived.   In  view  of  the  intended  study  of 
oscillations,  the  general  case  of  unsteady  motion  is  con- 
sidered.  Section  II  contains  the  study  of  the  equilibrium 
jet.   If  gravity  effects  are  neglected,  the  solution  can 
be  obtained  in  a  closed  form.   Corrections  due  to  gravity 
are  also  estimated.   In  section  III,  results  of  section  II 
are  discussed  and  compared  with  some  experimental  data. 
Section  IV  contains  a  study  of  the  unsteady  jet.   In  this 
case,  explicit  results  cannot  be  derived,  nor  are  they  de- 
sired due  to  the  incertitude  about  the  initial  jet  shape, 
but  a  satisfactory  approximate  treatment  is  obtained. 

As  to  the  equilibrium  jet,  our  results  agree  very  well 
with  the  available  experimental  data,  but  only  after  cor- 
rections due  to  gravity  are  taken  into  account.   See  fig.  I4., 
It  is  rather  surprising  that  gravity  ought  to  have  a  detect- 
able effect.   The  thin  jet  theory  leads  to  satisfactory 
results  only  if  the  mass  flux  exceeds  a  certain  critical 
value.   Below  this  critical  value,  the  mathematical  prob- 
lem for  the  equilibrium  jet  fails  to  have  a  solution, 
showing  that  in  this  case  the  thin  jet  assumption  does  not 
apply.   Concerning  the  jet  oscillations,  we  find  that  they 
will  be  damped.   The  frequency  of  these  oscillations  is 
found  to  be  of  the  order  of  2,000  cycles  per  second  when 
the  outside  pressure  is  normal.   If  the  pressure  is  re- 
duced this  frequency  may  drop  to  as  low  as  20-25  cycles  per 
second.   The  time  in  which  the  amplitude  of  oscillations 
will  be  reduced  to  half  its  original  value  is  found  to  be 
of  the  order  of  .03  second  and  is  independent  of  pressure. 

Other  Important  features  observed  in  connection  with 
conical  jets  are: 

1.  slow  period  oscillations  of  the  jet  envelope, 

2.  shattering  (or  dropletlzation)  when  the  mass  flux 
exceeds  a  certain  critical  value. 

The  theory  of  the  first  of  these  could  be  obtained  from 
the  general  equations  derived  in  section  I  provided  that 
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one  has  a  satisfactory  mechanism  for  the  production  of  a 
long  period  phenomenon.   It  has  been  suggested  that  air 
may  diffuse  through  the  envelope,  thus  resulting  in  a  dif- 
ferent volume-pressure  relation  for  the  air  enclosed  by  the 
jet  than  the  one  considered  in  section  IV.   This  would  re- 
sult in  a  different  formula  for  the  period,  but  the  same 
expression  for  the  half-time  of  such  oscillations.   As  to 
the  shattering,  the  techniques  employed  in  this  paper  are 
not  adequate  to  consider  this  phase.   A  study  of  it  is  now 
in  progress  and  will  be  discussed  in  a  future  report. 

The  author  has  been  unable  to  find  in  the  literature 
any  reference  to  unsteady  conical  jets.   As  to  the  steady 
case,  this  has  been  studied  by  T.  G.  Hodgkinson,  "Control 
by  Surface  Tension  of  a  Conical  Fluid  Sheet  Jet",  Porton 
Technical  Paper  No.  17^j  May  1950*   In  our  view,  Hodgkinson' s 
treatment  is  inconsistent  with  the  assumptions  he  makes. 

The  author  wishes  to  thank  Mr.  John  S.  Derr,  Jr.  for 
supplying  him  with  the  experimental  results  used  in  this 
paper. 
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I.   Formulation  of  the  problem, 


nozzle 


stream] 


Denote  by  r  ,  h  ,  a 
the  nozzle  mean  radius  at  orifice, 
nozzle  thickness,  and  nozzle 
angle,  respectively.   See  fig.  1. 
We  consider  the  flow  of  liquid 
emerging  from  the  nozzle  X'/ith 
velocity   v  .   If  p   is  the 
density  of  the  liquid  then 


(1) 


M 


2Trp  r  h  v 
'^  o  o  o 


is  the  mass  fliix  through  orifice. 

vie   assume  that  the  nozzle 
axis  is  vertical  so  that  even 
if  gravity  effects  are  considered, 
the  axial  symmetry  of  the  prob- 
lem is  preserved.   It  is  there- 
fore sufficient  to  confine  our 
attention  to  the  meridional 
section  of  the  jet.   Thus  we 
consider  the  liquid  contained 
between  two  axial  planes  making 
the  angle  A  (P  ,  see  Fig.  1. 

V/e  use  the  Lagrangian  des- 
cription of  the  liquid,  that  is, 
consider  particle  paths  which 
will  be  called  trajectories. 
What  one  actually  observes  is  not  the  trajectories  of 
individual  particles,  but  the  flow  pattern  at  a  fixed 
time,  which  may  be  described  in  terms  of  streamsurfaces, 
the  outer  streamsurfaces  being  the  boundaries  of  the  jet. 
We  shall  call  streamline  the  meridional  section  of  the 
median  streamsurface. 
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Under  the  thin  Het   assiimption  it  is  assuined  that  the 
velocity  field  is  constant  along  each  normal  to  the  stream- 
line.  It  then  follows  that  the  thin  jet  is  described  com- 
pletely by  the  streamline,  and  that  the  trajectories  of 
particles  leaving  the  orifice  at  the  same  time  are  identical. 
The  thickness  of  the  jet  measured  along  the  normal  to  the 
streamline  can  be  determined  from  the  incompressibillty 
assiimption. 

Denote  by  x(T:,t)  the  position  at  time   t   of  a 
particle  which  at  time  t  was  at  the  orifice.   Thus 

(2)  x(T;,t)  =  r(T:,t)a  +  z(T:,t)b 

where   a   and  b   are  unit  vectors  in  the  radial  and 
vertical  (downwards)  directions,  respectively,  and 

(3)  xCTjT)  =  r^a. 

For  a  fixed  'c,  x  is  defined  for   t^  <  t  <  T(t),  T(t)   being 
the  time  at  which  this  particle  reaches  the  axis  (i.e., 
r(T,  T(t))  =  0),  and  x   represents  the  trajectory  of  the 
particle  labelled  by  T.   Let   t  (t)   be  the  time  at  which 
the  particle  that  reaches  the  axis  at  time   t,  has  passed 
through  orifice.   Then,  for  t  (t)  <  t  <  t,  and  fixed  t, 
X   represents  the  totality  of  particles  that  are  observed 
at  time   t,  hence  yields  a  parametric  representation  of  the 
streamline,   t  (t)   is  Implicitly  defined  by  r(T;  (t),t)  =  0. 

Consider  the  mass   Am   of  fluid  that  flox>fed  through 
the  orifice  between  times  i     and  x  +  At, 

(I;)        ^"^  =  t^  '^^^'^- 

At  time  t  >  T  +  /^iT,  this  mass  occupies  the  position  on 
the  streamline  between  x(x  +  At,  t)   and  x(t,  t)   and 
has  an  area 


(5)       £is  =  rA^^  |x^|  At 
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in  common  with  the  median  streamsurface .   The  forces  acting 
on  J71  are : 

1.  force  of  inertia  =  -  Z\^4-t 

2.  force  of  gravity  =  ^mgb 

where   g  is  the  acceleration  of  gravity, 

3.  pressure  force  =  -/vp  £sq 

where   q   is  the  unit  normal  vector  to  the  streamsurface 
pointing  inwards,  and  £jp   is  the  excess  of  pressure  of  air 
enclosed  by  the  jet  over  the  outside  air  pressure, 

k-.      force  of  tension  =  2_^sTHq 
where  T  is  the  coefficient  of  surface  tension,  and  H  is 
the  sum  of  principal  curvatures  of  the  streamsurface.   The 
Nev/ton  equations  of  motion  for  the  mass  ija     now  yield 
on  using  (i|)  and  (5), 

(6)       14^^^  =  MgB'  +  2Trr|x^l(-^  +  2TH)q. 

We  have, 

(8)  H  =  .|x^r3(,^z^^  _  ,^r^^  +  r-lz^|x^|2)  '■'' 
and 

(9)  ZiP  =  P  -  P  » 

where  p  is  the  pressure  of  outside  air  (assumed  atmospheric) 

and  P   is  the  pressure  of  air  enclosed  by  the  jet  and  is 

a  function  of  the  volume   V  enclosed  by  the  jet.   If  the  jet 

it 

The  sign  of  H  is  so  chosen  that   H  is  positive  v;hen  the 

streamline  is  concave  toward  the  axis,  in  which  case  the 
force  of  tension  is  directed  inwards. 
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Is  either  vented,  or  does  not  close,  ^  =  0,   Observing 
that  . 

(10)  V  =  -rr  y    r^z  dT  =  V(t), 

T^(t) 

we  see  that  the  equations  of  motion  (6)  form  a  set  of  two 
inte^ro-partial -differential  equations .   These  equations 
determine  x  provided  that  we  specify  the  initial  jet, 
that  is,  x(t,0),  x.  ('t:,0),  and  the  conditions  at  the  orifice, 
The  latter  are  given  by  equation  (3)>  and  the  specification 
of  velocities  at  the  orifice  is 

(11)  v(t;,t)  =  X,  (TjT:)  =  v  sin  a  a  +  v  cos  a  b. 

T/  O  O 


II.   steady  flow  solution;  case  of  the  vented  ,iet . 

We  have  a  steady  flow  if  the  trajectories  of  particles 
leaving  the  orifice  at  different  times  are  identical.  That 
is,  for  all  u,  x,  and  t, 

x('t:  +  u,  t  +  u)  =  x('r,  t ) . 

Differentiating  with  respect  to  u,  and  setting  u  =  0  one 
gets 

x^  H-  x^  =  0 

and  this  shows  that  the  flow  is  steady  if  x   is  a  function 
of   t  -  1   only.   It  follows  then  that  the  trajectory  traces 
the  streamline. 

To  find  the  steady  flow  streamline  we  thus  assume  that 
x(t,  t)  =  x(t  -  t),  and  use  equations  (1-6)  in  which  x  ,  x 
can  be  now  replaced  by  -x.   and  x. .   respectively.   This 
equation  will  have  solutions  depending  on  t  -  t,  only  if 
qP   is  also  a  function  of   t  -  t.   Since   /j)   is  a  function 
of  t   and  not  of  t,  see  equation  (1-9)  and  (I-IO),  it 
follows  that  ^     must  be  constant,  which  we  assume. 
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Writing   t   in  place  of   t  -  T,  and  introducing  dimension- 
less  variables  defined  by 


(1) 


i  x(t)  =  r^x  (t  ) 


i  t  =  r.v:^t 


o"  o 


the  equations  (1-6)  now  become  (on  using  a  dot  to  denote 
differentiation  with  respect  to  t'   and  dropping  all 
primes). 


-3 


-1-  I  . ,2 


(2)    X  =  [J-b  +(2t,r  -  V  Ix  I  "-^rErz  -  zr  +  r~  z|x  |  ]  )  (za^  -  rb ) 


where 


[^   = 


(3) 


r^g 


i-iv 


\    '     -0 

[j.TTTr, 


V  = 


Mv 


The  initial  conditions  (1-3)  and  (I-ll)  become 
ik)  r(0)  =  1,   z(0)  =  0,   r(0)  =  sin  a,   z(0) 


cos  a, 


If  gravity  effects  are  neglected,  the  problem  can  be 
solved  in  closed  form.   If  Ap  ^  0,  the  solution  involves 
elliptic  functions  and  its  discussion  is  too  involved  to  be 
considered  here,   V/e  shall  therefore  consider  only  the 
case  of  a  vented  jet  (/vp  =  0,  i.e.  i,  =   0),  but  will  keep  the 
gravity  terms. 

Multiply  equations  (2)  scalarly  with  x'.   The  result  is 


X  •  X  =  -  (  |x|  )•  =  |J,Z  , 


(5 )       "^   "^  ~  2 

and  on  integrating  and  using  (i|.), 
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(6)  |x|^  =  r^  +  z^  =  1  +  2|iz. 

From  this  we  get  the  inequality   |z|  <  /I  +  2[iz    ,  and  as  a 
consequence,  the  inequality 

(7)  |z|  <  /l  +  2iJ,z  <  1  +  tit. 

This  shows  that  for  finite  t,  r,  r,  z,  and  z  remain  bounded. 
Now,  observing  that 

[(rz  -  zr)|xr^r  +  zlxT^jf  =  (zlxT^r)*  =  [— ^- .}' 

\  v/1  +  2m.  z  / 

the  b   component  of  equation  (2)  can  be  integrated  yielding 

(8)  (1 _^Z___)z  =  ^it  +  (1  -  v)  cos  a. 

/I  +  2!J,z 

The  remainder  of  our  considerations  depend  on  whether 
V  <  1   or  V  >  1,  and  we  consider  the  case   v  <  1   first. 
In  this  case  the  right  hand  side  of  (8)  remains  positive, 
and  since   z   is  bounded,  neither  of  the  factors  on  the  left 
hand  side  can  vanish.   Consequently,  they  preserve  their 
signs  at   t  =  0,  and  therefore 

(9)  z  >  0,   1 ^^      >  0. 

/I  +  2txz 

Incidentally,  using  (8)  and  inequality  (7)  we  now  get 

|a-t  +  (l-v)cosa=(l-    '^^  )z  <  ,/l  +  2iiz  -  vr 

/l  +  2m,z 

<  1  +  lat  -  vr, 
hence 

(10)  r  <  ^  [1  -  (1  -  v)  cos  a]  . 

If   V  =  1,  inequalities  (9)  still  hold  for   t  >  C,  hence  (10) 
holds  also.   ',/e  now  have   r  <  1,  and  since  this  contradicts 
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the  initial  conditions   r(0)  =  1,  r(0)  =  sin  a  >  0,  the 
problem  has  no  solution  for   v  =  1. 

Since   z  >  0,  z   can  be  chosen  as  the  independent  vari- 
able, and  we  consider   r  as  a  function  of  z.   The  equation 
that   r  satisfies  can  be  obtained  as  follows.   Vie  compute 
z   in  terms  of   r,  z   and   r  from  equation  (5)  and  substi- 
tute in  the   a  component  of  (2)  (which  does  not  involve   t 
explicitly).   On  solving  for   r,  we  thus  get   r  and  there- 
fore  z   in  terms  of   r,  z,  r,  z.   Now  differentiating  the 

relation  r  =  r'z   (primes  denoting  differentiation  i;ith 

ti 

respect  to   z)  we  obtain  r   in  terms  of   z,  r,  i,    hence  in 


terms  of 
using  z 


r,  z,  r,  z.   Finally,  r,  z   are  eliminated  by 


/  ^^  obtained  from  (6),  and 


r  =  r  z. 


The 


f  1  +  r 
result  of  these  manipulations  is 


(11) 


,"   =   «    M-r      +   vi/l    +    2M-Z 
vr 


1   - 


'2 

1   +  r  "^ 

1   +   2iaz 


/T 


ci\i.Z 


with  the  initial  conditions 


(12) 


r(0)  =1,   r  (0)  =  tan  a. 


Equation  (11)  can  be  integrated  explicitly  in  the 
case  when  ij.  =  0.  He   get  then 


r 
1  +  r 


and  on  integrating, 
(13)      /l   +  r'2 


cos  a 


1  -  vr 


1  -  V 


Equation  (13)  determines   |r  |   and  not   r   itself.   How- 

II 

ever,  in  view  of  (9),  r  <  0   for  all   z   (if   IJ.  =  0 )  and 
therefore   r   decreases  for  all   z.   With  this  information, 
the  integration  can  now  be  completed,  and  one  obtains 


aiK)   ri;?!'' 
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Uk)  r  =  -  -  c  cosh  (k  -  -) 

V  w 

where 

(15) 

Equation  (II4.)  shows  that   r  increases  to  a  maximum  r 

(16) 

attained  at 


^  (1-v)  cos  a    ^  ^  ^^3^-1  3^^  ^_ 
V      * 


r=--c=-[l-  (1-v)  cos  a] 


(17) 


z  =  z  =  ck 


and  then  decreases  to  zero  attained  at 


(18)      z=z=c(k+k') 


where 

(19)  k     =  cosh~      sec    p,      cos    p  =    (1-v)    cos    a. 


>z 


Figure  2 

The  behavior   of  r(z)   is  not  materially  changed  when 
tj,  >  0.   Indeed,  in  view  of  (9)   r   remains  negative  as 
long  as   r   is  positive.   Hence   r   decreases  to  zero  and 
then  become'S  negative.   Since  at   z  =  0,  r   has  now  a  lower 
value  than  in  the  case  [i  =   0 ,    the  curve  initially  runs 
below  the  one  obtained  before;  see  fig.  2.   As   r   decreases 
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II 
further,  r   will  eventually  get  closer  to  zero  (or  even 

change  its  sign)  resulting  in  a  flattening  of  the  curve. 

II 
However,  r   need  not  change  sign  before   r  becomes  zero 

(that  is,  before  the  jet  crosses  the  axis). 

In  case   v  >  1,  the  right  hand  term  of  equation  (8)  is 

initially  negative  but  changes  sign  at  the  critical  time 

^  (v-1)  cos  a    ^  remains  positive  for   t  >  t  .   Until 
on,  ^  o 

this  time  the  two  factors  on  the  left  hand  side  of  (8) 

preserve  their  signs  and  we  now  have 

vr 


(20) 


z  >  0, 


V/1  +  2|J,z 


-  1  >  0  for  t  <  t  . 

o 


Hence,  for  t<  t  ,  z   can  be  chosen  again  as  independent 
variable,  and  in  the  same  manner  one  arrives  at  equation 
(11)  for  r(z),  which  shows  that   r   remains  positive  and 
r   Increases.   If   tJ,  =  0,  one  obtains  thus  the  same  result 
for  r(z)   as  equation  (li|)  except  that  in  the  present  case 
c  <  0  and   r   increases  exponentially.   See  fig.  3« 


Figure  3 
If  the  gravity  term  is  included,  the  behaviour  of   r(z) 
changes  appreciably.   To  see  this,  we  observe  that  at  least 
one  of  the  factors  on  the  left  hand  side  of  equation  (8) 


must  vanish  at   t  =  t. 


Consider  one  of  these,  namely. 


vr  -  \/l   +  2nz  =  f(z).   For 
V  -  1  >  0,  f  (0)  =  V  tan  a 


t  =  0,  (i.e.  z  =  0),  f(0)  = 
■  \i,    which  we  assume  positive, 
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vhereas      f    (z)    =  vr      +  |J.'^(1   +   2txz)~-^^^   >   0     for     0  <   t   <   t 


L^n    +   2iiz)~^/^   > 


11  ,  -  o 

since      r      >  0.      Hence      f    (z)    >   0      for      0  <   t   <   t        so   that 

f(z)      increases,    hence  cannot    vanish.      We   conclude    then  that 

■  2  ' 

z(t^)    =  0,    and   since   then     r     =  1   +  2uz   >  0,      r        is   infinite 

O  It 

and  so  is  r  .   The  solution  cannot  be  continued  beyond  t  , 
for  then  z  would  have  to  decrease. 

Corrections  due  to  gravity  terms  can  be  estimated  in 
the  case   v  <  1.   Denote  by  r(z,  ^)      the  solution  of  (11) 
and  assume  that 

(21)  r(z,  ^l)  -^  r(z)  +  [xR(z) 

where   r(z)   is  the  zero  order  solution  given  by  (ll|-).   Sub- 
stitute in  (11)  and  expand  all  terms  in  power  series  in  \x, 
neglecting  terms  of  order  two  or  higher.   One  thus  obtains 
an  equation  for  R(z), 

(22)  c%'  +  2c  tanh  (k  -  -)r'  +  R  = 

-  [z  -  c  sinh  (k  -  -)  cosh  (k  -  -)] 

V  c  c 

with  the  initial  conditions 

R(0)  =  r' (0)  =  0. 

Since   sinh(k  -  — )   is  a  solution  of  the  associated  homo- 
geneous  equation,  (22)  can  be  solved  by  quadratures  and  the 
result  is 

(23)  R(z)  =  •^[z  -  c  sinh(k  -  ^)cosh(k  -  f )  +  2c  sec  a  sinh(k  -  ~) 

-  tan  a(z  sinh(k  -  — )  +  c  cosh(k  -  ~))^  > 

The  maximum  of   r(z,  ^i)   can  be  estimated  by  evaluating  it 
at   z  =  z.  We  get 

(2i;)      F(ti)  ■■  r(I)  +  txR(z)  =  1  -  c  -  ^  (tan  a  -  k). 


■HfelJ 
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The  point  where   r(z,  [i)      vanishes  can  be  estimated  by 
solving   r(z,  |i)  =  0  by  Newton's  method.   'Je  thus  get 

r  (z)  +  \xR    (2) 


or  explicitly 

(25)  z{[i)    =  elk  +  k'    + 


^  \  ,         (k+k' )cos(p-a)+3ec    p  Qin( p-a)-2sin   p' 

—  sin   p  cos    a  -    (k+k    )sin   a  -  -r^ 

[I  ^  1-v        _. 


In  a   typical   case,    a  =  30°,    v  =    .06      and     |i,  =    .00i|.      We   get 
r  =  3.100,      ^R(z)    =  -.025 

z'  =   16.12,      -IX      ,     .^  ^^^^    ,  =    .275. 

r    (2)    +  [iR    (z) 

Corrections  due  to  gravity  become  important  when  r  is  small. 
In  the  case  considered,  r(2z)  =  1,  whereas   ixR(2'z)  =  +.121, 
or  12°/o. 


III.   Comparison  with  experimental  results* 

Experiments  on  jets  produced  by  conical  nozzles  were 
performed  recently  by  J,  S.  Derr,  Jr.  at  Camp  Detrick, 
Frederick,  Maryland,  ^^i'ho  supplied  the  data  used  below.   The 
nozzle  used  had  the  following  characteristics: 

Radius  at  orifice,   r  =  .508 

'    o 

Nozzle  area  perpendicular  to  exit  section,  A  =  .O63I 

Nozzle  angle,  a  =  30 

Nozzle  thickness,  h  =  .0203. 

The  surface  tension  of  the  liquid  was  assumed  to  be   T  =  71, 

density,  C  =  1.   Three  flows  were  observed  v;ith  mass  flux 

M  =  21.5,  l8»'+j  1^1-. 2  respectively.   All  of  the  a- ove  data  are 

expressed  in  COS  units. 

M 
Since   v   =  7-7  ,  one  gets  for  v 


vj    ■■< 


J.ii  JX-iV  - 


irf^T::'!:, 
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(1) 


V  = 


liTrT  r  A 

0  0 

m2 


=  28.6  M 


-2 


using  the  above  data.   Hence  all  the  three  flows  considered 

fall  in  the  case   v  <  1.   If  one  assumes  that  the  initial 

jet  angle  is  equal  to  the  nozzle  angle,  the  sole  knowledge 

of  mass  flux  and  nozzle  characteristics  are  sufficient  to 

determine  the  theoretical  curves. 

Using  the  uncorrected  formulae  (III-II4.-I9)  we  now  get 

the  following  comparison  table,   (Lengths  are  expressed  in 

units  of  r  . ) 
o 


M 

21.5 

iS.ii 

li;.2 

V 

.0618 

.0851J. 

.li|l8 

r 

Computed 

3.05 

2*kk 

1.81 

Observed 

3.9i; 

3.50 

2Ak 

Error 

-25% 

-31V0 

-26% 

"z 

Computed 

7.23 

5.10 

2.88 

1 

Observed 

8.00 

7.39 

3.78 

Error 

-107o 

-327o 

-2kVo 

z 

Computed 

16,5 

11.70 

7.12 

Observed 

19.6 

15.87 

8.63 

Error 

-16% 

-26% 

-18% 

The  discrepancy  between  the  data  thus  computed  and  the 
observed  data  is  too  big  to  be  accounted  for  by  corrections 
due  to  the  gravity  terms.   The  cause  of  this  discrepancy 
lies  in  that  there  must  be  some  sticking  of  the  liquid  to 


.  i  !  L  l 
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the  nozzle  causing  the  initial  jet  angle  to  be  different 
from  that  of  the  nozzle.   Indeed,  a  further  scrutiny  of 
the  experimental  data  shows  that  this  angle  is  larger  than 
the  nozzle  angle  and  increases  when  the  mass  flux  decreases, 
and, as  formulae  (III  -  II4.-I9)  show,  r,  z,  z  are  sensitive 
to  the  variation  of  this  angle. 

To  determine  a  theoretical  curve,  a,  v,  and  n  must  be 
known.   These  are  related  by  the  condition 


(2) 


v^ 
M- 


2T 


f  Sr^h^  cos  a 


cos  a 


obtained  from  (I-l)  and  (II-3),  involving  only  static 
quantities.   Here  we  denoted  by  a     the  nozzle  angle 
(equal  to  30°  in  the  present  case)  which  now  is  assumed 
different  from  the  initial  jet  angle   a.   Since   a  cannot 
be  measured  accurately,  it  is  preferred  to  determine   v,  [i, 
and   a  by  requiring  that  predicted  values  of   r  and   z 
agree  with  the  experimental  values.   Using  the  equations 
(II-2I4.)  and  (11-25)  for   r  and   z,  we  thus  obtain  together 
with  (2)  three  equations  for   v,  [i,    and   a.   For  the  three 
flows  considered  these  yield 


Case  (M) 

a° 

V 

\i' 

1  (21.5) 

31^.3 

.0559 

.ooi^iB 

2  (iQ.k) 

36.3 

.0720 

.00553 

3  (li^.2) 

39.8 

.1386 

.01117 

r(z)   computed  with  these  data  for   a,  v,  \i     are 
graphed  on  fig.  k-,    showing  an  excellent  agreement  with  the 
experimental  data.   In  case  2  and  3  the  discrepancy  never 
exceeds  2°/o  of  measured  data.   In  case  1  there  is,  at  some 
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points,  a  divergence  of  about  10°/o,  but  this  is  undoubted- 
ly due  to  a  poor  extrapolation  of  measured  data  which  was 
necessary  to  obtain  the  value  of   z. 

As  to  the  case   v  >  1,  no  experimental  data  are  avail- 

1    2    T 
able.   As  formula  (II-3)  shows,  v  >  1  when  •^  t.v   *^  vT"  ' 

o 

which  would  occur  when  v  <  83.6  cm/sec   or  M  <  5*28  g/sec. 

It  is  doubtful  whether  in  this  case  the  thin  jet  assumption 

leads  to  coherent  results. 

IV.   Theory  of  oscillations  of  unsteady  jets. 

If  the  initial  jet  conditions,  x(t,0)  and  x^(t,0), 
(see  section  I),  deviate  from  the  equilibrixom  jet  discussed 
^n  section  II,  an  oscillatory  motion  is  expected.   The 
mechanism  leading  to  such  a  motion  was  already  discussed  in 
the  introduction. 

In  general  it  is  very  difficult  to  solve  the  initial 
value  problem  for  the  system  of  equations  (1-6).   Fortunately 
such  a  solution  is  not  even  desired,  for  it  is  not  natural 
to  prescribe  specific  initial  conditions  beyond  the  fact 
that  they  deviate  from  the  equilibrium  configuration.   An 
approximate  theory  is  obtained  from  the  consideration  of  the 
energy  equation  obtained  by  dotting  equation  (1-6)  with  v 
and  integrating  over  the  jet  envelope  (i.e.,  with  respect 
to  1)    at  a  fixed  time.   The  interval  of  integration  is 
(t^(t),  t). 

We  consider  various  terms  separately.   Denoting  by   v 

and   V   normal  and  tangential  components  of  velocity,  we 
s 

have 

t  ft 

(1)   M   /    X.  ,  -  V  dt  =  ^  M  /    (v^^  +  v^)dT 

T:(t)   ^*  '^o 

o  o 

=  I  M  /   v^^d^  .  yi   M  ;   v/dT)  .  I  M(v2^  -^  -   v2^). 


■L  i. 
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Here,  v    is  the  normal  component  of  the  initial  jet  velocity 
and  V  ^   is  the  normal  component  of  the  velocity  at  the  tip 
of  the  jet.   Next,  it  is  easily  verified  that 

(2)       2Tr    /    r|x JhC  .  li^dT  =  -  ^fl^ 
T^Ct) 

where   S   is  the  area  of  the  streamsurface, 
S  =  2tt    /    rlx  I  dT, 


T^(t) 


and  that 


where   V(t)   is  the  volume  enclosed  by  the  jet,  see 
equation  (I-IO).   If  we  denote  by  G(v)   the  internal  energy 
of  the  air  enclosed  by  the  jet, 

ih)  G(V)  =  /  P(7)d-T   , 

V 

the  two  last  terms  in  (1-6)  contribute  on  integration 
-  1^  (+  G  +  pV  +  2TS). 

Combining  the  above  V7e  thus  get  on  integrating  (1-6), 

.t 

(5)       ft  ^1  ^    ■  ^n^^'^  ^   ^^^■^  ^  P^  ^  2TS] 


'T.(t) 

dt  (t) 


°  t 


T^(t) 


=  0. 


The  first  term  in  the  bracket  can  be  viewed  as  the  kinetic 
energy  of  the  normal  motion,  whereas  the  others  -  as  the 
potential  energy. 
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If  the  jet  were  in  equilibrium,  the  normal  component  of 
the  fluid  velocity  would  vanish  and  the  tangential  com- 
ponent would  satisfy 


1    2 

according  to  equation  (II-6)  (written  in  original  units). 
Furthermore,  since  v  would  be  practically  constant  and 
equal  to  v  ,  the  equation  of  the  jet  could  be  written  in 
the  form 

(6)  r  =  f(z),   f(0)  =  r^, 
with  the  parametric  representation 

(7)  z  =  z(T),   r  =  f(z('r)),   r^^  +  z^  =  vf   . 

Since  our  object  is  primarily  the  study  of  deviations  from 
the  equilibrium  jet,  it  is  plausible  to  ass\ame  that  at  each 
time  the  jet  is  approximately  in  equilibrium,  so  that  it  has 
a  representation  of  the  form  (6)  and  (7),  now  depending  on 
time  through  some  parameter   R(t),   That  is 

(6')      r  =  f(z,R),   f(0,R)  =  r   =  a  constant  in  time. 

In  this  approximation  the  last  term  in  (5)  vanishes. 
The  normal  velocity  component  is 

f^R 

(8)  v^  =   .-.^-  -2- 

and  from  (6')  one  gets   v   =0.   Furthermore,  if  the  para- 
metric representation  (7)  is  so  chosen  that   z(0)  =  0,  the 
time  that  a  particle  takes  to  flow  through  the  jet,  i.e., 
t  -  T-Ct),  is  equal  to  the  value  of   t  for  which  f(z,R) 
vanishes.   This  value  depends  on  R  and  will  be  denoted  by 
t(R).   Hence, 

dT  (t)    ,  , 

(9)  °^   =  It  ^^  ■  '^^^^^  =  1  -  t:  (R)R. 


J  .r-...- 


■;:  V .'.:'.  y  . 


•  r^-     .:j;>iu 


'I     » 
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Likewise,  the  final  value  of   z,  i.e.   zCtCR))   will  depend 
on   R.   Using  the  notation 

T(R)   fp^ 

(10)  M  /      2 — ^  dT  =  m(R), 

'0     1  +  f'*^ 

(11)  G(v)  +  pV  +  2TS  =  q(R), 

,      f^(2(T(R)),  R)  v^„  2 

(12)  iM  ii— ^ =  1  M(-^)   =  s(R), 

"^^    1  +  f ''^(z(T(R)),R)    ^     R 

equation  (5)  now  becomes  an  equation  for  R(t),  namely 

^■^^^      It  ^1  '^^^  +  q)  +  s(l  -  T*R)R^  =  0. 

Equation  (13)  for  R(t)   actually  describes  the 
response  of  a  thin  jet  which  at  each  time  is  constrained 
to  take  the  form  prescribed  by  (6'),  and  which  is  traversed 
by  an  incompressible  fluid  with  velocity  v  .   The  des- 
cription will  be  satisfactory  if  equation  (13)  has  an  equil- 
ibrium solution  R  =  R  =  constant,  and  for  R  ,  f(z,R  ) 

o  '  o'     '  o 

is  the  equilibrium  jet  found  in  section  II,  whereas  for  R 
differing  slightly  from  R^,  f(z,R)   differs  little  from  the 
equilibrium  shape. 

The  manifold  of  solutions  of  (13)  can  be  discussed 
without  reference  to  the  specific  choice  of   f(z,R),  ex- 
cept for  some  mild  assumptions  on  continuity  and  differ- 
entiability of  functions  appearing  in  (13)-   We  need  not  be 
concerned  with  these  here,  but  shall  assume,  in  accordance 
with  further  applications  that  there  exists  a  unique  value 

R   for  which 
o 


(111) 


q'(R^)  =  0, 


I  iimq"  -  s^  >  0. 


Carrying  out  the  differentiation  in  (13)  we  get,  on  dividing 
by   R, 


:;'.:vWj 


#J        -. 


tf^'i 
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(15) 


mR 


R  +  sR  +  (2^  -  ST' )R^  +  q'  =  0. 


It  is  Immediately  verified  that   R  =  R   is  a  solution  of 

''  o 

(15)  and  we  shall  now  show  that  any  solution  of  (15)  v;hich 
gets  close  to   R   performs  damped  oscillations  about   R  , 
In  a   t_  interval  where   R(t)   is  monotonic,  a  unique 
inverse  t(R)   is  defined.   Introducing  u(R)   defined  by 


(16) 


u(R)  =  R(t(R)) 


one  gets  from  (15)  a  first  order  equation  for  u. 


(17) 


u'  =  - 


+  su  +  ( 


m' 


ST' )u' 


mu 


The  point   R=R,u=0   is  a  singular  point  of  this 
equation,  and  an  easy  verification  shows  that  it  is,  pro- 
vided that  assumption  (lU-)    holds,  a  vortex  point.   Con- 
sequently, the  solutions  of  (1?)  near   R  =  R  ,  u  =  0   spiral 
about  this  point,  see  fig.  5»   The  reinterpretation  of  this 
"phase  space"  diagram  in  the  "physical  space"  leads  to 
damped  oscillatory  motion. 


locus  u'=0 


u'<0 


u'<0 


u'>0 


Fip;ure  5 

Since  R  satisfies  a  non-linear  equation,  its  "period" 
(time  elapsed  between  two  successive  maxima)  and  its  "half 
time"  depend  on  initial  data.   One  obtains  an  estimate  of 
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the  period  and  the  half  time  by  assuming  a  motion  of  small 
amplitude  and  linearizing  the  differential  equation.  Set- 
ting R  =  R   +  eR,  and  linearizing  (15)  one  now  gets 

(18)  m(R^)R  +  s(R^)R  +  q"(RQ)R  =  0. 

Equation  (18)  has  solutions  of  the  form  exp(-K  +  /^i)t, 
with  period 

(19)  T  =  4^  =      ^^  ..    , 


and  half-time 


(20)     e  =  i  =  f 

(in  the  above,  m,  s,  and  q"   are  evaluated  at   R  =  R  ), 

The  above  considerations  indicate  that  the  jet  will 
perform  a  damped  oscillatory  motion,  provided  that  the  im- 
portant assumption  (1[|.)  is  satisfied.   Actual  estimates 
of  T   and   0  depend  on  the  choice  of  f(z,R).   Since  our 
theory  is  only  approximate  there  is  no  criterion  for  the 
best  choice  of  f .   Undoubtedly  a  satisfactory  result  would 
be  obtained  by  choosing  for  f   the  function 

r    (1  -  R)r  cos  a  _, 

(21)      r  =  f(z,R)  .=  ^ g cosh[cosh"   sec  a 

Rz 1 

(1  -  R)rQ  cos  a-" 

which  for  R  =  v   reduces  to  the  equilibrium  solution  for 
the  jet  envelope  derived  in  section  II  (eq.  II-II4.).   This 
function,  which  is  obtained  from  {ll-ll\.)    by  replacing   v 
by   R  has  also  the  advantage  that  its  slope  at   z  =  0   is 
tan  a,  (no  matter  what   R   is),  so  that  it  represents  at 
each  instant  an  equilibrium  jet  corresponding  to  some  mass 


•t    >;> 
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flow,  emerging  with  the  same  angle. 

To  use  (21)  for   f(z,R)   would  lead,  however,  to  com- 
plicated calculations.   It  is  felt  that  not  much  is  lost  by 
using  a  "jet"  shape  which  is  conveniently  simple,  but  either 
has  an  equilibrium  area  equal  to  that  of  a  jet  actually 
considered,  or  encloses  the  same  volume  in  equilibrium. 
For  the  sake  of  simplicity  we  therefore  assume  that  the  jet 
is  a  sphere  of  radius   R(t)   attached  to  the  nozzle.   See 
Fig.  6. 
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^jet 
envelope 


Its   equation   is 


(22) 


r   =   f(z,R)    =     i/r^   -    (z   -    C)^    ,    where   C   = 


■M  -  rj. 


The  desired  parametric  representation  is 

v 


and 


i  z  =  C  -  R  sin[sin"""^  ^  -  "5"  '^J 


-1  C 

R  '     R 


^  V 

I  r  =  R  costs m   ^  -  -^  tJ , 


[  T(R)  =f  (f-^  sln-l|), 
/         o 


z('r(R))  =  R  +  C. 
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Since 
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V  =  Ik  '"  C  '"'^'  -  ""^  ' 


R+C 


S.  =d_  (2^  {—   f  ^f.2  dz)  =  2Tr(RlC) 


dR 


2 
C   ~  ' 


one  obtains 


q'(R)  =  -PV  +  pV  +  2TS'  =  ^^^^„'^^^' 


-p(v)  +  P  +  ^i 


using  the  definition  (11).   The  equilibrivun  radius   R   is 
therefore  determined  as  the  solution  of  equation 

P(V)  =  P  +  ^  . 


Assuming  the  adiabatic  relation  for  the  air  enclosed  by 
the  jet,  P  =  kV~',  we  now  get  for  the  equilibrium  value  of 


q  » 


,2- 


U3)   ,",R^,  =2M^;xll,p.!f,  -kTl 

R  J 


_  ttR(R+C)' 

c 


C[R(R+C)+r^]       ^     R^^l 


evaluated  for  R  =  R  .   Using  the  definitions  (10)  and  (12) 
for  m(R  )   and   s(R  )   we  find 


(21+) 
(25) 


r3r. 


m  = 


MR  


+  (1  +  ^)(§  +  sin"^  g)  , 


2C 


2'  '2 


R' 


s    2  ^  C  ^  • 


Even  the  simple  model  {2.2)    leads  to  rather  cumbersome 
formulae  for  the  period  and  the  half-time.   A  further  simp- 
lification is  obtained  on  expanding  all  formulae  in  powers 
of   r^/R.   We  get 
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(23')  q"(R^)    =   12TrYpR(l    +  ^    (^)^   +    ...) 

-  l6Tr(3Y-l)T(l   -i^^f^   if)      -    ...), 
(2V)  m=|^   (1   .i   (^)2.    ...), 

(2^')  s    =   2M(1    +  ^    (r^)      +    ...). 

In  the   application  to  be   considered,    -^  <  -i-   .      Hence, 
omitting  the   terms    involving     r        will  result   in   an   error 


o 


of  less  than  0.5  /o  in  case  of   q",  and  of  about  5  /o  in  the 
evaluation  of  m  and   s.   Incidentally,  setting  r  =  0 
amounts  to  taking  in  place  of  the  "jet"  described  by  (22), 
a  full  sphere  touching  the  nozzle  at  one  point.   Setting 
r^  =  0  and  substituting  (23'}-(25')  in  (19)  and  (20),  we 
get  the  approximate  formulae 

(27)      «  =  ^  • 

o 

As  an  application  we  consider  the  three  cases  dis- 
cussed in  section  III.   As  was  already  indicated,  we  could 
require  that  either  the  area  of  the  sphere  is  equal  to  the 
area  of  the  actual  jet,  or  the  volume  of  the  sphere  is  equal 
to  the  volume  enclosed  by  the  actual  jet.   Each  of  these 
criteria  will  lead  to  a  slightly  different  value  for  the 

equilibrium  radius.   V/e  use   R   for  the  radius  deter- 

o  _ 

mined  by  the  first  of  these  and  R   for  the  radius  deter- 
mined  by  the  second.   One  gets,  in  COS  units. 
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Case 

M 

^0 

V 

S 

% 

«o 

1 

21.5 

3k.l 

68.5 

92.)+ 

2.13 

2.51+ 

2 

18. i; 

292 

la.o 

61+. 6 

2.29 

2.11+ 

3 

11+. 2 

225 

12.1+ 

26.7 

1.50 

1.1+3 

Using  R^,  in  (26)  and  (27),  we  finally  get 


e 

T 

Case 

p=l   atmos. 

p=0 

1 

.035 

.1+5  X  10"^ 

.053 

2 

.031^ 

.1+5  X  io"3 

.01+8 

3 

.030 

.1+5  X  10"^ 

.039 

T/iThen  p  >  .001  atmos,,  the  last  two  terms  in  the  bracket 
of  (26)  are  unimportant,  and  the  formula  for  the  period  re- 
duces to 


(28) 


1.9 


1/  p 


where   A   is  the  area  of  the  exit  section  of  the  nozzle. 
Ivhen  the  pressure  is  below  this  value,  formula  (2  9)  becomes 
inapplicable,  but  as  the  table  shows,  we  get  about  20  to  25 
cycles  per  second.   The  half-time,  on  the  other  hand,  is 
independent  of  pressure,  and  is  of  the  order  of  .03  sec. 
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